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Abstract 

We discuss neutrinoless double beta decay and lepton flavor violating decays such as fi — > cy in the colored 
seesaw scenario. In this mechanism, neutrino masses are generated at one-loop via the exchange of TeV-scale 
fermionic and scalar color octets. The same particles mediate lepton number and flavor violating processes. 
We show that within this framework a dominant color octet contribution to neutrinoless double beta decay 
is possible without being in conflict with constraints from lepton flavor violating processes. We furthermore 
compare the "direct" color octet contribution to neutrinoless double beta decay with the "indirect" contribution, 
namely the usual standard light Majorana neutrino exchange. For degenerate color octet fermionic states both 
contributions are proportional to the usual effective mass, while for non-degenerate octet fermions this feature 
is not present. Depending on the model parameters, either of the contributions can be dominant. 
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1 Introduction 



The question regarding the origin of neutrino mass remains one of the most pressing problems in particle physics [T] . 
In particular, knowledge of whether neutrinos are Dirac or Majorana particles is of paramount importance. If 
neutrinos were to be Majorana particles, lepton number would be violated. Since lepton number is an accidental 
symmetry of the standard model of particle physics, one would like to know the scale of lepton number violation 
(LNV), if any, leading to the generation of neutrino masses. In addition, the peculiar pattern of neutrino masses and 
mixing demands an explanation which requires the extension of the standard model, and one would like to probe 
the existence of new particles usually considered within such extensions. While neutrino oscillations and beta decay 
experiments continue to improve our knowledge of the low energy neutrino mass matrix, we have to look at data 
from a variety of complementary avenues such as neutrinoless double beta decay (0v[3j3) experiments [2j|5], lepton 
flavor violation searches (6j[7j, as well as collider experiments to augment our understanding of the mechanism 
responsible for neutrino mass generation. On the experimental side, there are realistic prospects for order-of- 
magnitude improvements in the search for the 0v/3j3 half-life (81. Moreover, there has been an order-of-magnitude 
improvement on the bound for fi — > ej recently (6] . 

While one usually assumes that the exchange of light Majorana neutrinos is the leading contribution to neutri- 
noless double beta decay [9] , it is well known that a plethora of alternative intermediate scale theories exists (where 
the intermediate scale is TeV, few hundred GeV and even few tens of GeV), which not only predict neutrinoless 
double beta decay, but can also saturate the current bounds on the process, or can lead to sizable rates in current 



and future experiments 10-22 . Some of the interesting alternative approaches are based on seesaw models 19 - 21 



left-right symmetric theories 18 , i?-parity violating supcrsymmetry 11 ■ 15 and so on. See the recent reviews 23 



and 24 for summaries of the particle physics aspects and the experimental situation of 0vj3/3, respectively. The 
fact that heavy TeV-scale particles can lead to the same contribution to 0v/3f3 as sub-eV scale neutrinos is easily 
understood from looking at the particle physics amplitude in the standard interpretation of light neutrino exchange: 



A, 



t—' (p 2 ) — n 



(2.7 TeV) ~ 



(1) 



Here we have inserted the current limit of about 0.5 eV on the effective neutrino mass M ee — | m i\i with U 

being the lepton mixing matrix, rrii the neutrino masses, and (p 2 ) ~ 0.01 GeV 2 3> m 2 for the relevant momentum 
scale of the process in the neutrino propagator. Thus, if heavy particles with mass much larger than (p 2 ) are 
exchanged in the process, for instance two scalars or vectors with mass Mi and one fermion with mass M%, they 



will contribute with 
TeV regime. 



l/(Mf M 2 ) to the amplitude, and can saturate the current limit if their masses are in the 



In this paper we work with a particular variant of the seesaw mechanism 25 which has TeV-scale scalar and 



fermionic color octets. The fermionic octets have a Majorana mass term leading to lepton number violation in the 
theory. Consequently, one expects that they will directly generate neutrinoless double beta decay. Neutrino masses 
are forbidden at tree level, however they are generated at one-loop level, via the exchange of color octet fermions 
and scalars. This model has been dubbed the "colored seesaw" model 26 . Since the new particles introduced 
in this model are color octets, sizable cross-sections at the LHC can be expected, connected with the spectacular 
feature of lepton number violation 27 . In addition, the lepton flavor violation associated with the gauge invariant 
Yukawa couplings of the octets 



2.s 



implies that not only the neutrino sector is non-trivial in what regards flavor, 
but also the charged lepton sector is. Hence, decays like fx — > provide additional tests of the scenario 28 , 
and the branching ratios depend on the same set of parameters as the neutrino mass matrix. Here, we discuss 
neutrinoless double beta decay and lepton flavor violation mediated by these color octet states, and show that 
they, being TeV-scale particles, can give a large and in fact saturating contribution in both processes. Note that 
the current collider limits on their masses are 1.92 TeV for the octet scalar (311, and about 1 TeV for the octet 



fermions 32 



The reason why these contributions have not been considered in the literature so far is that usually one assumes 
minimal flavor violation (MFV) to avoid large flavor changing neutral current (FCNC) effects [33]. Then, the 
coupling of the charged member in the weak doublet scalar octet to a quark q is proportional to m q /v, with v being 
the vacuum expectation value (VEV) of the standard model Higgs. Therefore, the amplitude of 0^/3/3 receives a 
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Figure 1: Neutrino mass generated at the one- loop level in the colored seesaw scenario. 



suppression factor of mr u d /v and is completely negligible 26 . We depart from the assumption of minimal flavor 
violation in this work, and hence a sizable rate of 0^/3/3 can come from the color octet states. In this framework, we 
study the interplay of neutrino mass and mixing, lepton flavor violation and neutrinoless double beta decay, and 
encounter several interesting and general features, which are characteristic for the correlations between the various 
phenomenological sectors. 

One particularly interesting feature is the interplay of direct and indirect contributions of the color octets 
to Ou/3/3. With "direct contribution" we mean the new contribution with exchange of heavy octets, while with 
"indirect contribution" we mean the usually considered light Majorana neutrino exchange mechanism, where here, 
however, the neutrino masses are generated by the octets via loops. While in general these two contributions 
depend differently on the parameters, we identify situations in which they are both proportional to the effective 
mass, namely when the fermionic octets are degenerate in mass. In other cases, they depend differently on the 
flavor parameters. In all cases, the neutrino exchange mechanism can be dominating or sub-leading, depending on 
the masses of the octets and the quartic coupling governing the interaction between the color octet scalar and the 
standard model Higgs boson. 

The outline of the paper is the following: we start with the basics of the colored seesaw scenario in Section [2j 
Following this, in Section |3j we discuss the contribution of color octet states in neutrinoless double beta decay and 
lepton flavor violating processes in general, and consider in some detail the most simple case of two degenerate 
color octet fermions in Section |4j The general case of three non-degenerate color octet fermions is the subject of 
Section [5] before we conclude and summarize in Section [6] 



2 Colored Seesaw Scenario 

In the colored seesaw mechanism [26|[27| , the particle content of the standard model is extended with a color octet 
scalar $ and color octet fermions ^j, which transform as 

*~(8,2,+l) ; *i~(8,l,0). (2) 

While $ is charged under SU(2)l and U(1)y, the fields are singlets under these gauge groups. Note that while 
the model just contains one additional scalar, we need at least two additional fermions in order to generate at least 
two massive light neutrinos which are required to explain current neutrino data. We will consider the cases i = 1, 2 
and i = 1, 2, 3 in this paper. The relevant Lagrangian corresponding to the new sector is 

-C v = YfL a ia 2 + iM^Tr(^) + A $ff Tr($t H) 2 + h.c, (3) 

where H is the standard model doublet Higgs and L a is the lepton doublet of flavor a = e, /it, r, i.e., greek indices 
correspond to flavor and roman indices to physical mass states. We have considered A$jy and M^ i as real. Since 
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the scalar $ is a color octet and SU(3) C symmetry must remain unbroken, it must have a zero vacuum expectation 
value. Hence, it is evident from the above that at tree level one cannot generate the neutrino mass operator 
0( LH ^ H ) [29] . However, neutrino masses will be generated at one-loop level through the mediation of colored 
octet scalars and fermions, as shown in 26,27] and illustrated in Fig. [I] The one- loop neutrino mass matrix is 

\ „,2 ^<&H va i v l3i 



M' 



16tt 2 



(4) 



where the loop function X is given by 



Xi =X(M^,Mi Si ) = M* t - 



Ml+M* lnfe) 



{Ml-Mlf 



(5) 



Here, is the mass of the scalar octet, and v = 174 GeV is the doublet Higgs VEV. The neutrino mass matrix in 
Eq. @ can be diagonalized by the unitary 3x3 Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix U, 
where according to our convention U^M V U* = M d = diag(mi, mg, m^), with Wj being the light neutrino masses. 
Using this and Eq. Q, one can express the Yukawa coupling matrix Y v in terms of light neutrino masses, mixings, 
color octet masses and a Casas-Ibarra 30 matrix 1Z as: 



'16tt 2 1 



U JM$KJ(I d )- 1 



where I d = diag^i,^,^). In general, we can express the complex orthogonal matrix 1Z as 



(6) 



C2C3 

TZ = I -C1S3 - S1S2C3 
S1S3 - C1S2C3 



C2S3 s 2 
C163 - S1S2S3 S1C2 
-S1C3 - C1S2S3 cic 2 



(7) 



with S{ = sin#;, c$ = cos^ (i = 1,2,3), where 9±, #2 and O3 are arbitrary complex angles. Since the color octet 
scalar field <E> is charged under SU(3) C as well as SU(2)l and U(1)y, it can interact with the quark fields of the 
standard model. Note that the scalar field $ can be decomposed as 



$ = A 



$0 



(8) 



where the X A , A = 1, ... ,8, are the Gell-Mann matrices. The Lagrangian corresponding to the interaction between 
the quarks and the scalar $ is 



Co = d R K D &Q L + u r kuQ l § + h.c. 



(9) 



where Ku,d are the corresponding Yukawa couplings. Going from the flavor basis to the physical basis of the quark 
fields (Ul.r and P>l,r are the bases rotation matrices for the up and down quarks, respectively), the Lagrangian is 
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We choose the basis where the up-quark mass matrix is diagonal, i.e., the up-quark mixing matrices are Ul,r. 
Using this, the above equation simplifies to 
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Figure 2: Neutrinoless double beta decay me- 
diated by light Majorana neutrinos. 
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Figure 3: Neutrinoless double beta decay me- 
diated by color octet scalars and fermions. 



After describing the basic Lagrangian in this section, we now discuss the different processes where these color 
octet states can participate. In particular, we concentrate on the novel color octet contribution to neutrinoless 
double beta decay and on lepton flavor violating processes like \i — > &y. We will show in particular that the color 
octet framework for neutrino mass generation can easily saturate both the fi — ¥ ej limit as well as the limit on 
neutrinoless double beta decay. 



3 Neutrinoless Double Beta Decay and Flavor Violation 



The standard light neutrino contribution to 0i//3/3 is given in Fig. [2j It is proportional to the effective mass 



Mee = £E£ 



We note in this paper that the scalar and fcrmionic states will also contribute in this lepton 



number violating process, as shown in Fig. [3] The conceptually interesting part is that there is a direct and an 
indirect contribution from the octets. The direct one is the new diagram from Fig. [3j while the indirect one is the 
neutrino exchange mechanism in Fig. [2] It is called indirect because the light neutrino masses are generated by the 
octets in the loop diagram given in Fig. [T] The relevant vertices for the direct contribution of the octets to Oj//3/3 



A a/3 U 



Pl {kuDl) 1x $id and r i/ ei e L ($ + )** l , 



(12) 



so that the relevant effective operator for Qv(5j3 is (uuee\C cB " (x)\dd), where C eS " (x) is 
{Y » l)2 MlM^ Xaf}XlS ( Pr ^d d r)ii ~ Pl(kuD l ) u ) d p ] [u 7 (Pr(4,D r ) u - Pl^uD l ) u ) d s ] (e L e c L ) . (13) 

For simplicity, and for illustration, we consider the case where Ku -C Kd and concentrate on the right-chiral part 
only. We denote the contribution coming from the quark states as y\ x — yfi/3, where yn = (k^D^h, and the 
factor p can come from the hadronization procedure. The inverse half-life of 0vj3f3 is given by the well-known 
expression 

™ P vli \- (y?y 



1 

1/2 



= G 0v \M v r) v + M^rj^\ , where r) v = 



and ?7$* = 



(14) 



Here M. v and A^$* are the nuclear matrix elements corresponding to the light neutrino exchange (indirect con- 
tribution) and to the exchange of color octets (direct contribution). We will focus here for definiteness on the 
neutrinoless double beta decay of 76 Ge, for which a half-life limit of Ti/ 2 > 1.9 x 10 25 yr exists [2j. The phase space 
factor is Gqu = 7.93 x 10 _15 yr _1 . Alternatively, one can give the above expression as: 
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Figure 4: Lepton flavor violating process l a — > Ipj in the colored seesaw scenario. 



In the above, JCqu = GQ V {M.^m p ) 2 , P 2 — m e m p ^"" and Gf = 1.166 x 10 5 GeV 2 is the Fermi constant. The 
parameter P 2 describes the difference between the long- and short-range contributions. 

It is also instructive to compare the amplitudes for neutrinoless double beta decay on the particle physics level: 



Here ^4 is the amplitude for the direct contribution (color octet exchange) and Ai the indirect contribution due to 
the light Majorana neutrino exchange. We deal with a nuclear physics problem here, so that a typical size of the 
momentum is (p 2 ) « (100 MeV) 2 . Comparing different amplitudes on the particle physics level is very often an 
excellent approximation to compare their relative contribution to Qvpp. With the help of the complex orthogonal 
matrix 1Z one can re-write the color octet amplitude in terms of the PMNS matrix: 

One can immediately see from this expression that for degenerate fermion octets, i.e., for M^ k = Mqr the elements 
of 1Z drop out of this expression and the amplitude is proportional to M ee = J2 m i^ei- Therefore, if the fermion 
octets are degenerate, both the direct as well as the indirect contributions to 0^/3/3 are directly proportional to 
the effective mass M ee . As a result, if due to any cancellation the contribution from light neutrino exchange goes 
to zero, the corresponding contribution from the octet exchange also vanishes identically. However barring these 
cancellations due the Majorana phases, we will see that the relative importance of the two contributions may widely 
vary depending on the model parameters. 

Let us point out that in some sense Qvf3fi and the generation of neutrino masses can be decoupled in this model. 
It can be seen from Eqs. ^ and ^ that two independent sources of lepton number violation exist. The first term 
of Eq. |3) may be used to assign lepton number to the color octet particles, and the last two terms then break it. 
If A<j>£r is zero, we can assign lepton number either to $ or to ^i. In the first case (assigning lepton number to $) 
the quark couplings to $ in Eq. |9| provide LNV; in the second case (assigning lepton number to the colored 
fermion mass term gives LNV. If now = 0, the one-loop neutrino masses in Eq. Q vanisfQ but still there are 
non-vanishing contributions to 0vf3(3. 

Besides neutrinoless double beta decay, the color octets can actively participate in different lepton flavor violating 
processes, like fi ej, r — > and r — > /17, as shown in Fig. The branching ratio corresponding to l a — > Ipj 



is 26 28 



(17) 



1 Since lepton number is not conserved, higher order diagrams will lead to very small neutrino masses. 
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where the function T{xi) has the following form, 



l-6xi + Zxj + 2x1 - Qx i Mxi) 
I2(xi - l) 4 



(18) 



with Xi — M^./Mq and a em being the fine structure constant. As for the case of the Ov/3/3 amplitude, one can 
write the branching ratio of \x — > in terms of the neutrino parameters and the matrix 1Z as 



Br(/i — > ey) 
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(19) 



For the simple choice of TZu — 8™, the above expression reduces to 



Br(/z — ¥ &y) 



3a c 



(16^ 2 ) 2 



(20) 



In fact, one can check that the branching ratio is independent of 7Z if the fermion octets are degenerate and 1Z 
is real. Since 71 embodies our lack of knowledge of the seesaw-scale physics which cannot be determined from 
low-scale data, we can immediately see that for such a class of model parameters the theory is fully determined 
leading to unambiguous prediction for lepton flavor violating decays. Currently, the best limit on lepton flavor 
violation comes from the MEG collaboration, which gives [6] 



Br(/i ej) < 2.4 x 10 



-12 



(21) 



at 90% C.L. In what follows we will discuss some simple examples for the interplay of neutrino mixing, neutrinoless 
double beta decay and lepton flavor violation. While this is not a complete analysis, some very interesting and 
general features arise. 



4 Neutrinoless Double Beta Decay and Lepton Flavor Violation with 
Two Color Octet Fermions 

In this section, we discuss neutrinoless double beta decay and lepton flavor violating processes, considering the 
minimal case with two degenerate color octet fermions. Thus the matrix 1Z depends on only one complex parameter 
and can be taken as 

/ o o \ /yr^ 2 - u \ 

ft(for NH) = VI-lo 2 -u) and 7£(for IH) = I uj ^Jl-u 2 . (22) 
\ u VT^j 2 J \ / 

In our analysis we consider a real u and hence, — 1 < to < +1. Note that since we need one heavy fermion for each 
light neutrino, this results in the lightest neutrino mass being zero. We first discuss the normal hierarchy and then 
the inverted hierarchy scenario. 



4.1 Normal Hierarchy 

For the normal hierarchy scenario and two color octet fermions with degenerate masses, i.e., Mq, i = (and hence 
I, = T) for all i, the Yukawas can be expressed as 



Y v = W-^-C/diag(0,V^2:,\/^)^diag(Vl rT ,v / Z- 1 ). (23) 

V *<S>H V 
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Figure 5: Half-life T X j 2 of 76 Ge against the mass of the color octet fermion M$. The mass of scalar octet has 
been fixed as M$ = 1.96 TeV. In the left plot, the red, blue and orange lines correspond to A$# = 0.15 x 10~ 8 , 
\<s>h = 0.1 x 10~ 8 , and A$# = 0.7 x 10 -9 , respectively. The yellow region is experimentally excluded. The black 
line corresponds to the Heidelberg-Moscow limit. The brown line represents the light neutrino contribution to 
ncutrinolcss double beta decay divided by 50 2 , for the neutrino oscillation parameters mentioned in the text. In the 
right plot the red, blue and orange lines correspond to the three different nuclear matrix elements — 600.38, 

M.<t,^ = 377.59, and M.^ = 188.79, respectively. The coupling has been set to \$> H — 10 -9 . Both figures are for 
normal hierarchy. 



Expressed in terms of the elements of the PMNS mixing matrix U, the light neutrino mass eigenvalues rrii and the 
quartic coupling A$#, the corresponding half-life for 0vf3/3 can be written as 

' -KoJ±+ „, !L 2 .1^ S \rn 2 Ul 2 + mM\ (24) 



As discussed in the previous section, since we have considered degenerate fermions, the expression for the half-life of 
neutrinoless double beta decay is independent of the parameter w, although the Yukawa couplings depend strongly 
on it. We also stress again that for the case of degenerate octet fermions both the light neutrino contribution and 
the color octet contributions share the same proportionality factor |M ee | = \m2U 2 2 + TnaU%\. Hence, if due to any 
cancellation between neutrino parameters the light neutrino contribution becomes zero, the same cancellation will 
set the color octet contribution to zero, too. However, for M ee ^ the direct contribution from the color octets 
can be taken independently of the neutrino mass scale. This is because there is an additional parameter A$# in 
addition to the Yukawa coupling involved, and hence one can suitably adjust the two to get small enough neutrino 
masses required for normal hierarchy and yet get a very large contribution to 0^/3/3. In particular, by choosing the 
coupling A$# to smaller values, the color octet contribution can dominate over the light neutrino contribution. We 
have shown the relative comparison between the light neutrino and color octet contributions in Fig. [5j where we 
have plotted the variation of the half-life T\/ 2 vs. the mass of the color octet fermion My. The details of the figure 
are as follows: 



We used the neutrino oscillation parameters Am\ x = 7.5 x 10 5 eV 2 , Am§ x = 2.3 x 10 3 eV 2 , 6x2 — 33° 



23 = 42° and 6*13 = 8°, see Ref. 35 for a recent global fit. For this choice of parameters, the effective mass 



is M ee = 0.003 eV. The elements U e 2 and U e 3 have been considered real. 

The mass of the scalar octet has been fixed as M$ = 1.96 TeV, so that it satisfies the present bound on the 



mass of color octet scalars 31 . However, in 31 the bound has been derived assuming negligible coupling 
with quarks, i.e., production in the gluon fusion channel. In our case, the color octet scalar has interaction 
with the quarks, hence the mass bound may be weakened. 

The black line represents the Heidelberg-Moscow limit [2] for neutrinoless double beta decay and corresponds 
to T1/2 = 1.9 x 10 25 yr. The yellow region is thus excluded, while the blue region is the allowed one. 
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Figure 6: Branching ratio of /i — > ej as a function of the octet fermion mass for normal hierarchy. The mass of 
scalar octet has been fixed as M$ = 1.96 TeV. The red, blue and orange lines correspond to three different values 
\$h = 0.15 x 10~ 8 , A$# = 0.1 x 10~ 8 , and A$# = 0.7 x 10 -9 . The yellow region is experimentally excluded. The 
black line corresponds to the present bound from the MEG experiment. 

• The brown line represents the standard light neutrino contribution to O/v/3/3 scaled by a factor 50 2 , i.e., 

K-Ov (l^) x 552 > wnere we have taken P 2 ~ (251.49) 2 MeV 2 . The factor K-ov depends on the nuclear 
matrix element and can be obtained using ICq^ = GQ^{M^m p ) 2 . We have taken y\ x = 1 and adopted the 



nuclear matrix elements from 13-15 , where M. v = 2.8, and following Eq. (A. 8) in 15 , A4$* = 377.59, if 
the pion exchange and usual one and two nucleon mode are considered. For this value of M$^, ICoi, — 9.95 x 
10~ 10 yr" 1 GeV 2 and the standard contribution is K.^ (^1f )~ 2 = 3.38 x 10 29 yr. These O^/3/3-matrix element 
values are the ones which have been evaluated for the case of short range i?-parity violating SUSY diagrams, 
which is closest to the setup considered by us. In particular, the operator structure which corresponds to 
A' m diagrams is similar to our case. Hence the information regarding nuclear matrix element for i?-parity 
violation can be utilized in our scenario. 

The red, blue and orange lines correspond to the values A$# = 0.15 x 10~ 8 , A$# = 0.1 x 10 -8 and A$# = 
0.7 x 10~ 9 , respectively. The half-life decreases with A$#. 

• In Fig. [5j we have also shown the variation of Ti/ 2 with the mass of the color octet fermion, considering the 
variation of nuclear matrix element in the interval 188.79 — 600.38. 

The figure clearly shows that the direct octet contribution can dominate the contribution from light neutrino 
exchange easily, which, we reiterate, has been shown by dividing by a factor of (50) 2 . It is also clear that the direct 
octet contribution can easily saturate the current limit on the half-life of 0^/3/3, even for normal hierarchy. Note 
that the prediction for 0vf3f3 with light Majorana neutrino exchange for normal hierarchy is very low compared to 
that for inverted hierarchy, and much below the reach of the next generation Qvj3(3 experiments. However, since 
in the colored seesaw model one can have very large 0z//3/3 even for normal hierarchy, 0v/3f3 can no longer be used 
distinguish between the neutrino mass hierarchies. 

Having discussed lepton number violation, we turn to lepton flavor violation now. We are interested in exploring 
if the color-octet fermions can simultaneously produce a 0v/3(3 rate large enough to saturate the current limit, and 
at the same time saturate the experimental bound for fi — > ej. The branching ratio for this process has been given 



in Eq. ( 19 ). We stick to real uj and degenerate fermions and hence the branching ratio becomes independent of the 



parameter uj. It has the following form: 

B* - n) - s^Ffjj |^ (^)' \rn^ 2 U; 2 + m3 U^ . ,25, 

Like in neutrinoless double beta decay, the color octet fermions can also give a significant contribution to this 
process. The result has been shown in Fig. [6j where we have given the variation of the branching ratio against the 
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mass of the color octet fermions, considering three different values of A&h ■ The summary of the figure is as follows: 



• The neutrino oscillation parameters as well as the mass of scalar octet have been set to the previously 
mentioned values, already used in the 0vf3f3 study. For 12 = 33°, 6> 23 = 42°, 13 = 8°, Ato^ = 7.5 x 1(T 5 eV 2 
and Amf^ = 2.3 x 10~ 3 eV 2 , we have (rn^U^V^ +fn^U e3 U^) = 7.1 x 10~ 12 GeV. The elements of the mixing 
matrix have been considered real. 

• The black line corresponds to the present bound coming from /i — ¥ ej searches at MEG experiment, i.e., 
Br(> ->■ &fj = 2.4 x 10~ 12 . 

• The red, blue and orange lines correspond to the values \<$>h = 0.15 X 10~ 8 , = 0.1 x 10~ 8 , and A$# = 
0.7 x 10 -9 , respectively. 

It is clear from the figure that one can rather easily saturate the current lepton flavor violation limits. Note that 
in principle the expression \m 2 U e2 U* 2 + fn 3 U e ^U* 3 \, and thus the branching ratio for /i — > cy, can vanish if the 
mixing parameters and CP phase conspire. This would also be possible when we considered complex 1Z. 

If the color octet contribution saturates both the Heidelberg-Moscow bound, as well as the present bound on 
[i — > ej decay coming from MEG experiment, then the mass of the scalar and fermionic octets, and the coupling 
\<£H should satisfy simultaneously the following two equations: 

A$h = 1.63 x 10 4 yu 1 j 



Ml 

We have represented the above two conditions in Fig. [7j where we have shown the variation of A$# against the 
mass of the color octet fermions M*. As before, the mass of the scalar octet and /Coi/ have been set to be the same 
value as for Fig. [5] For simplicity, we have considered the elements of the PMNS mixing matrices to be real. In 
Fig.[7J the gray region is excluded both from 0vf3f3 and fi ej. The red and blue lines in the left panel correspond 
to Bi(fi — > ej) = 2.4 x 10 -12 , and the 0v(3f3 saturating bound T\/2 — 1.9 x 10 25 yr, where the Yukawas between 
scalar octets and quarks have been considered 0(1). The same exercise has been repeated for the right panel of 
Fig. [7] with a different y\ x factor. Note that with the decrease of y\ x , one will require an additional suppression 
in \m 2 U e2 U* 2 + m 3U e 3U* 3 \ in order to simultaneously saturate the MEG and Heidelberg-Moscow limit. This can 
come from cancellation between the phases. Also, interestingly Fig. [7] indicates that if only lepton flavor violation 
is saturating, the color octet fermions can be within the reach of LHC. However, inclusion of a saturating Oj//3/3 
demands the color-octet fermion mass to be higher. 

4.2 Inverted Hierarchy 

We discuss briefly the relative comparison of neutrinoless double beta decay and /1 — > ej when the light neutrino 
states follow the inverted hierarchy. The half-life for Ov/3/3 can then be expressed as 

1 v ( 1 £11 1671-2 Vl ty2 , 7T 2 I 2 , 07 ^ 

= )Co» p. + jrr- tjt^o -or- -if m 1 U el +m 2 U e2 , (27) 



T 1/2 W \P 2 M 9 M*G%v*\* H lJ 
and the branching ratio for fi — > ej is given by 

Br(, - - jj^^ {^fj + m.tW&f . (28) 

Since, in general even for light neutrino states, the contribution for inverted hierarchy is larger than the contribu- 
tion for normal hierarchy, the same feature holds for the colored seesaw scenario with degenerate fermions. One 
can obtain a significantly large neutrinoless double beta decay contribution with relatively small values of quark 
Yukawas, and/or larger color octet scalar and fermion masses. For completeness, we show two figures to support 
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M^(GeV) 



Figure 7: The coupling against the mass of the color octet fermions My. The mass of scalar octet has been 
fixed as M$ = 1.96 TeV. The gray band is experimentally excluded by fi — > e 7 and Oz//3/3. Left panel: the red line 
represents the present limit obtained by the MEG experiment for n — > e 7 . The blue line represents the Heidelberg- 
Moscow bound. The factor y\ x has been considered to be 1. Right panel: The blue and purple lines correspond to 
the Heidelberg-Moscow bound, while the y\ x factor has been considered as 0.4 and 0.2 respectively. The red and 

pink lines correspond to the MEG limit, while we have considered a suppression factor 0.08 and 0.0197 from the 

1 1 2 

factor \m 2 U e2 U* 2 + m 3 U e3 U* 3 \ . The figure has been generated considering normal hierarchy. 



this feature. See Figs. [8] and [9] for more details. All parameters except the ones explicitly mentioned in the captions 
are the same as for the previous figures. 

We would like to end this section by discussing the particular relations which the different lepton flavor violating 
processes share among themselves in the special case considered here. The ratio of the branching ratios for different 
lepton flavor violating processes are 



Br(r -> e 7 ) = (U e2 U; 2 m 2 + U e3 U* 3 m 3 ) Br(r -> /17) = (U^U^m-i + U&U^ma) 

Br(/x -> e 7 ) (U e2 U* 2 m 2 + U e3 U* 3 m 3 ) ' Br(^i -> erf) (U e2 U* 2 m 2 + U e3 U* 3 m 3 ) 

for normal hierarchy, and 

Br(r -> e 7 ) = (U^U^mi + U e2 U* 2 m 2 ) Br(r -> M7 ) = {U^U^rm + U^U; 2 m 2 ) 

Bi(fi^ej) {U el U; imi + U e2 U; 2 m 2 ) ' Br(/i -> erf) (U el U* imi + U e2 U* 2 m 2 ) 



(29) 



(30) 



for inverted hierarchy. Note that, since for the two octet fermion case with real 1Z there is no additional parameter in 
the theory apart from the ones measured in low energy experiments, one can get an exact prediction for these ratios 
of branching ratios in terms of oscillation parameters. Hence, experimental measurement of some or all of them 
can be used to check if the lepton flavor violation is solely due to the neutrino mass generation mechanism or not. 
The present bounds are Br(r — > jjrf) < 4.5 x 10~ 8 and Br(r — > erf) < 1.2 x 10~ 7 , respectively [7|. Disregarding the 



possibility of cancellations in the numerator or denominator, the above expressions (29) and (|30|) are expected to be 
of order one. Since experimentally the limits obey Br(/i — > e 7 ) <C Br(r — > (i<y) as well as Br(/i — > e 7 ) <C Br(r — > e 7 ), 
the limits on t — > /i 7 and r — > e 7 are therefore automatically obeyed if the limit on fi — > e 7 is obeyed. 



5 Neutrinoless Double Beta Decay with Three Color Octet Fermions 

After discussing lepton number and lepton flavor violating processes for two degenerate color octet fermionic 
states in the last section, we now turn to the three generation scenario, dropping in addition the assumption of 



degenerate octet fermions. It is evident from Eq. ( 16 ) and the discussion following it that in this case the color 
octet contribution and the light neutrino contribution of 0vj3f3 do not share the same proportionality factor M ee 
anymore. This particular feature brings up the possibility that, even if the light neutrino contribution becomes 
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Figure 8: Half-life of Ov/3/3 in 76 Ge and the branching ratio for ^ — > ey as a function of the octet fermion mass M$ 
for inverted hierarchy. The mass of scalar octet has been fixed as M$ = 1.96 TeV. The red and blue lines correspond 
to \<j,H = 0.15 x 10~ 8 and A$ir = 0.1 x 10 -8 , respectively. In the left panel the light neutrino contribution has 
been shown in gray, without any scaling factor. While for the left plot y\ x = 0.05, the right panel is generated 
without any additional suppression factor. The blue areas are allowed, the yellow ones forbidden. 



zero due to cancellation between the terms with mi, m,2 and m-3, the color octet contribution can be non-zero, and 
even significantly large. We explicitly show this specific feature for one case. 

In general, the color octet contribution will have a significant dependence on the phases of the matrix 1Z. We 
do not address this issue in the present work, as we will encounter extreme cases even when IZtj = 8ij, in which 
case the particle physics amplitude of Qvfifi is given by 

_ 16tt 2 ( x ^ m,U^ \ 

We will stick to this simple case throughout the rest of this section, as very interesting features arise already at 
this stage. 

In what regards neutrino mixing, we try to keep things as simple as possible, and study a somewhat minimal 
deviation from tri-bimaximal mixing. Denoting sin #13 = A, the PMNS mixing matrix is now 



/ j_ x -is 



PMNS 



7S - 75 7S + 75 e " -73 I eM ' ^ +5) ) ■ (32) 

V6 V6 \/3 %/6 s/% 



Note that we allow for a complex PMNS matrix in what follows. In this case, the amplitude for the light neutrino 
contribution to 01/(3/3 is 



For the amplitude of 0^/3/3 mediated by the color octet fermions and scalars, we have 



(33) 



. Vu 16tt 2 / 2mi m 2 e l2a m 3 AV 2 ^ . 



The branching ratio for the process /1 — > ey is given by 

Br(/U — > ey) oc 
where the proportionality factor is | 11 



(35) 



3 AL» 4 G% Mi ' 
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Figure 9: The coupling X^h against the mass of the color octet fermions My for inverted hierarchy. The gray 
band is experimentally excluded by /j, — > and 0v(3/3. The mass of scalar octet has been fixed as A/$ = 1.96 TeV. 
The red and blue lines represent the saturating contribution in MEG experiment and Heidelberg-Moscow bound. 
The factor y\ x has been set to 0.05. 



Focussing on the ratio between A and Ai, and therefore on the relative size of the direct and indirect contri- 



butions to Ov/3/3, see Fig. 10 We have plotted there for the normal and inverted mass ordering the ratio of the 
amplitudes as a function of the smallest neutrino mass for different values of . It is obvious that even for this 
simple example that the ratio of the amplitude can be very large or very small, corresponding to the dominance of 
one of the contributions. The details of the figure are as follows: 

• The mass of the color octet scalar M$ has been set to 2 TeV, while the masses of the color octet fermions 
My. have been varied inside the interval [0.9, 1.1] TeV. No particular form of hierarchy has been considered 
between the color octet fermions. The random variation inside this interval mainly assures that there is no 
exact degeneracy between the three fermions. 

• The Yukawa coupling yn has been varied inside the interval [0.001, 1.0]. Additionally, all phases in the PMNS 
matrix have been varied in the interval [0, 2ir]. 

• The solar and atmo sphe ric mass-squared differences, as well as #13 have been varied inside their presently 



allowed 3er intervals 



35 



. The typical momentum scale has been been set to (p ) ~ (100) MeV 



• The differently colored regions in this figure correspond to different \$h values as shown, all satisfying the 
MEG limit Br(> -> erf) < 2.4 x 10~ 12 . 

• Note that as increases, the ratio A/Ai decreases. The large increase in A/Ai for low values of m\ in 
normal hierarchy is an artifact of phase cancellation in Ai . 

Another illustrative way to visualize the different contributions is to define an "effective mass" for the direct 



octet contribution. Noting that from the indirect amplitude the usual effective mass M ee 36 37 is obtaine d b y 



multiplying A\ with (p 2 )/G 2 F , we can define the "color effective mass" as %J-A, where A is given in Eq. (34) 



We can plot now both the standard effective mass M ee and its analogous expression ^jQ-A as a function of the 
lightest neutrino mass eigenvalue. The two plots are given in Fig. [IT] We see that the usual phenomenology can 
be significantly modified. For instance, in the inverted hierarchy (negligible 7713) one expects in the standard case 
M ee <; 0.05 eV. The direct contribution from the octets does approach 1 eV, and hence (for the simple example 
considered here), can be used to cut in the parameter space of couplings and masses. Note also that the predicted 
Oisf3(3 is very large even for normal hierarchy and almost comparable to that for inverted hierarchy, as pointed out 
earlier for the two octet case. 
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Figure 10: The ratio of particle physics amplitudes for normal hierarch y (l eft panel) and inverted hierarchy (right 



panel) as a function of the lightest neutrino mass for non-zero sin 2 # 13 34 . The parameters used are given in the 
text. The colored areas give the allowed regions using the current MEG limit for /i — > ej. Recall that the standard 
effective mass, or Ai, is observable in the next few years for values of mi <; 0.3 eV. The black and red lines give the 
upper or lower values of the correspondingly colored areas. Note that the areas for different values of overlap. 
The lower red and yellow lines are below the scale of the axis. 



6 Final Remarks and Conclusions 

There are two issues regarding our framework which we now address. The first one deals with dropping the MFV 
hypothesis. Strong constraints on FCNC processes exist, for example from scalar color octet exchange in K°-K° 



mixing, from b — > or the electric dipole moment of the neutron, see e.g. Ref. 33 . We note that the octet 
contribution to neutrinoless double decay that we consider in this work depends only on the coupling of the scalar 
octet $ with an up- and a down-quark. One can convince oneself that in all possible FCNC diagrams this coupling 
never appears on its own. For instance, in K°-K° mixing diagrams or in b — > sj it appears together with couplings 
involving 2nd and 3rd generation quarks. Constraints coming from the electric dipole moment of the neutron can 
be avoided by setting a possible phase to small values (in analogy to the SUSY CP problem). While this is not a 
completely satisfying situation, we nevertheless note that in the limit of only the coupling to up- and a down-quarks 
being non-zero, we face no phcnomcnological problem. In addition, neutrinoless double decay is the only place in 
which that coupling appears on its own and hence it is the only place where it can directly be constrained. 

Another point is the strong hierarchy between and the Yukawa coupling of up- and down quarks. Radiative 
corrections might spoil this hierarchy, for instance one might have diagrams in which the quartic HH coupling 
is mediated by quark loops. However, in the same limit as above, namely only the coupling of $ to up- and a 
down-quarks being non-zero, this diagram is suppressed heavily by (m Ui d/v) 2 and causes no problem. 

To sum up, in this work we have discussed neutrinoless double beta decay and lepton flavor violation in fi — > 
for the so-called colored seesaw scenario. In this model, the color octet scalars and fermions generate Majorana 
masses for the light neutrinos via one-loop diagrams. Since these states have non-trivial charge under SU(3) C and 
the electroweak gauge group, the same set of fields can directly participate in neutrinoless double beta decay and 
lepton flavor violating processes. 

Studying only simple examples, we have already found interesting features: it is for instance possible that the 
octet states saturate the limits on both (j, —¥ ej and neutrinoless double beta decay. If the octet fermions are 
degenerate in mass, then the contributions to Ov/3/3 from the octets and the light neutrinos are both proportional 
to the effective mass M ee , their relative importance depending on the model parameters. 

It is conceptually interesting that the color octets imply a direct and an indirect contribution to neutrinoless 
double beta decay: the direct contribution is the one considered here for the first time, namely the short-range 
exchange of octet scalars and fermions. The indirect contribution is the standard, long-range one with the exchange 
of light Majorana neutrinos. These neutrinos are generated at loop-level by the octets and it is interesting to 
compare those two contributions. Extreme cases are easily possible, in the sense that both contributions can be 
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Figure 11: Left: the usual plot [36] of the effective mass against the smallest neutrino mass. Right: the "color 
effective mass" %s L A 1 which has to be compared with the effective mass in the standard diagram for 0v(3(3, as 

F 

a function of the lightest neutrino mass. For the right plot, = 10~ 8 . The normal ordering is given in blue 
and the inverted ordering in yellow. The darker areas are valid when the oscillation parameters are fixed to their 
best-fit values (for better visibility only best-fit values for the octet contribution are given), the brighter areas for 
the 3tr ranges, and sin #13 has been fixed to 0.013. In the right plot, the lower yellow line gives the minimum value 
for inverted hierarchy. Recall that the standard effective mass M ee is observable in the next few years for values 
of nil ~ 0.3 eV, and the current limit on the half- life corresponds to about 1 eV for M ee , which roughly is also the 
limit for 

either dominant or negligible. 
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